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ABSTRACT 

Two-dimensional thermal models for automated tape placement (ATP) of thermoplastic 
composites neglect the diffusive heat transport that occurs between the newly placed tape and the 
cool substrate beside it. Such lateral transport can cool the tape edges prematurely and weaken 
the bond. The three-dimensional, steady state, thermal transport equation is solved by the 
Green’s function method for a tape of finite width being placed on an infinitely wide substrate. 
The isotherm for the glass transition temperature on the weld interface is used to determine the 
distance inward from the tape edge that is prematurely cooled, called the cooling incursion A a . 
For the Langley ATP robot, A a = 0.4 mm for a unidirectional lay-up of PEEK/carbon fiber 
composite, and Aa = 1.2 mm for an isotropic lay-up. A formula for A a is developed and 
applied to a wide range of operating conditions. A surprise finding is that A a need not 
decrease as the Peclet number Pe becomes very large, where Pe is the dimensionless ratio of 
inertial to diffusive heat transport. Conformable rollers that increase the consolidation length 
would also increase A a , unless other changes are made, such as proportionally increasing the 
material speed. To compensate for premature edge cooling, the thermal input could be extended 
past the tape edges by the amount A a . This method should help achieve uniform weld strength 
and crystallinity across the width of the tape. 
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1. INTRODUCTION 

The strength and high-temperature properties of carbon-fiber thermoplastic composites make 
them ideal for many structural applications. However, fabrication with these composites can be 
labor-intensive and costly. Automated tape placement (ATP) is a process designed to reduce this 
cost. In ATP, thermoplastic composite tape is placed on the previous layer and welded by heat 
and pressure to build up a laminated structure, as shown in Fig. 1. 

Mathematical models have been developed to control and optimize the ATP process. Most of 
these models (refs. 1-10) are two-dimensional, as shown, for example, in Fig. 2. They restrict 
the heat transport to the x-y plane, and neglect the heat transport in the z-direction. In refs. 5-10, 
this two-dimensional approximation is justified partly by stating that the Peclet number Pe » 1 , 
where Pe is the dimensionless ratio of inertial to diffusive heat transport. The present 
modeling study should determine if sufficiently large values of Pe do, indeed, justify the two- 
dimensional approximation. 

If the applied heat q 0 is uniform across the tape width and zero outside, as shown in Fig. 1, a 
strong thermal gradient occurs in the substrate along the side edges of the tape. This thermal 
gradient can prematurely cool the tape edges and weaken the bond by causing heat to flow from 
the hot tape to the cool substrate beside it. 

The main purpose of this study is to help achieve uniform weld strength and crystallinity across 
the width of the tape. For this purpose, the model will determine the prematurely cooled distance 
inward from the tape edge (the cooling incursion A a ) over a large range of processing 
parameters. This range includes the present operating condition of the ATP robot at Langley 
Research Center. The approach is to use the Green’s function technique to obtain an analytic 
solution of the three-dimensional thermal transport problem shown in Fig. 1. Related studies at 
the leading edge of a new tape and at the trailing edge of a dropped tape are presented in refs. 8 
and 9. 


2. THREE-DIMENSIONAL MODEL 

Figure 1 shows the steady state configuration where the laydown/consolidation head is fixed and 
the material moves with speed U beneath it. The Cartesian coordinate frame is fixed with the 
head, and its origin lies at the center of the nip line. The thermal input q Q is directed at the nip 
line and is constant over the width of the tape (-a < z < a) and zero elsewhere. Far away from 
the nip, the supporting tool and the head are both maintained at temperature T v , which is 
usually room temperature. The thermal input q Q is constrained so that the temperature on the 
weld interface always decreases passively to the glass transition temperature T g at the 
downstream end of the consolidation head x h , where the pressure is released, as shown in Fig. 

2. This constraint is enforced in the two-dimensional limit as a — > 

T(x h ,0,z) = T g 


a — > °o 


( 1 ) 



3 


where z remains finite. In spite of this constraint, three-dimensional effects will cool the weld 
interface near the tape edge to T g in a shorter distance than x h . 

The supporting tool and the head are taken to have the same thermal conductivity as the average 
thermal conductivity of the composite substrate. This assumption enables the configuration 
shown in Fig. 1 to be approximated by a line thermal source of length 2 a that is fixed along the 
z -axis in an infinite moving conductor. The three-dimensional thermal transport equation then 
becomes 


P C p U %~ K n^- K 22 ^T- K 33 = ^ 0 <5(x)<5(y)[5(s + a)- S{z - a)] (2) 

where 

p density, kg m" 3 

C p specific heat at constant pressure, J kg 1 K 1 

U material speed in the positive x -direction, m s' 1 

T temperature, K 

K n , K 22 , K 33 tensor conductivities along the x -, y -, and z -axes, W m 1 K 1 
q 0 thermal input power density at nip, W m 1 

8 (x), 8 (y) Dirac delta functions, m" 1 
S(z) Heaviside unit step function 

a Half width of tape, m 

Dimensionless (primed) quantities are defined as follows: 


T'= 


T-T^ 

T^T~ 


(3a) 




2 7c(K u Kj' 





(3b) 

(3c) 

(3d) 

(3e) 
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The partial derivatives become 


and (2) can be written 


d_ 

dx 

d_ 

dy 

d_ 

dz 


1 d 


x h dx' 




1/2 


K, 


( K ' 


1/2 


V^3, 


a 

ay 

_a 

dz 


where 


dr d 2 T' d 2 T' 
dx' dx' 2 dy' 2 


d 2 T' 

w 


= q'(x',y',z') 


(3f) 

(3g) 

(3h) 


(4a) 


q'{x',y',z')=27rq Q '8(x')S(y')[s(z'+a')-S{z'-a')] (4b) 

subject to the boundary condition that 7”— > 0 as x ,y' , or z' — » °° • In dimensionless variables, 
the constraint (1) becomes 


lim 

a '— > oo 


r(i,o,z') = i 


(5) 


where z' remains finite. In order to solve (4) in the infinite domain, an adjoint Green’s function 
H' (x', y', z' I x 0 ' , y 0 ' , z 0 ') is defined that satisfies the adjoint equation 


dH' d 2 H' d 2 H' 
dx' dx' 2 dy' 2 


d 2 H' 


8 {x'-x 0 ')8{y'-y 0 ')8 {z'-z 0 ') 


(6) 


where H' also vanishes at infinity. Multiplying (4a) by H' and (6) by T' , subtracting, and 
integrating by parts over the infinite x',y',z' -domain gives 


T'{) 


o > 


Jo .Z< 


')= J dx J dy' \dz'q'{x',y',z')H'{x',y',z' \x 0 '. 


Jo 


(7) 


Equation (6) can be solved by the method of Fourier transforms to obtain the adjoint Green’s 
function 
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H'(x',y',z'\ x 0 


expl 


-V) 2 + (y-y 0 ') 2 + (z'-V) 2 . 

1/2 


J) 

4; r[ 

(x'-Xq'Y + {y'-y 0 'Y + iz'-z 0 'Y_ 

1/2 


( 8 ) 


If the x',y',z' -variables are interchanged with the x 0 ',y 0 ',z 0 ' -variables, equations (7), (4b), and 
(8) become 


'{x',y',z')= J dx 0 ' j dy 0 ' ^dz 0 'q'{x 0 ',y 0 ',z 0 ')H'{x 0 ',y 0 ',z 0 '\x',y',z') (9a) 


where 


/' (*0 ' ’ ^0 ' . Zo ') = 2n do ' S (-*0 ' ) s ( (z 0 '+«' ') ■ - S (zo ~ a ')] 


(9b) 


H’{x 0 ',y Q ',z 0 '\x',y',z')-- 


expl 

- b{[ {x 0 '-x'f + (y 0 '-y'Y + (z 0 '-z') 2 . 

1/2 

+ (.r 0 '-r') 

f) 

4 n\ 

(*o -v') 2 + (y 0 -y'Y + (z 0 '-z'Y . 

1/2 


(9c) 


Substituting (9b) and (9c) into (9a) and integrating gives 


T'{x\y\z') = ^-]dz 


expl 

r b \ 

[[x' 2 +y' 2 + (z 0 '-z') 2 . 

1 , 2 - 4 ) 

[ 

i2 , ?2 . ( i »\2 

X +y +U 0 -z) \ 

1/2 


( 10 ) 


As mentioned, the dimensionless thermal input q 0 ' is determined by the constraint (5). If (10) is 
substituted into (5) and the identification is made that 



( 11 ) 


where u' is a dummy variable of integration and 


K„(b) is a hyperbolic Bessel function, then 


do 


exp (b)K 0 (b) 


( 12 ) 


The dimensionless temperature field (10) becomes finally 
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r(x',y',z') = 


1 1 



[[x a +y 2 +{z,'-zf 

] ,, 2 - 4 ) 

2 exp {b)K 0 {b)\ 

| r 

x' 2 +y' 2 + (z 0 '-z') 2 ] 

1/2 


(13) 


Note that this solution depends on only two dimensionless parameters: b and the dimensionless 
half-width a . 


3. ISOTHERMS ON THE WELD INTERFACE 

As shown in Figs. 1 and 2, the weld interface occurs on the plane y'= 0 . Equation (13) can be 
used to plot the isotherms of dimensionless temperature 7" on this plane. (Computations are 
actually done on the plane y'= 0.001 to avoid the infinite temperature along the nip line 

(x'= 0, y'= 0, - a'< z'< a').) 


The interface isotherms are plotted in Fig. 3 for a = 1.5 and three values of b . The value 
b = 10 corresponds to the present operating condition of the Langley ATP robot. Cooling along 
the edge of the tape (z'=a'=1.5) is clearly evident. In the absence of this cooling, the isotherm 
r=l would lie on x'=l across the full width of the tape, in accordance with the constraint 

(5). 


The T'= 1 isotherm, which corresponds to the glass transition temperature, is used to quantify 
the dimensionless cooling incursion A a ' , as follows: The dimensionless distance from the tape 
edge (z'=a'=1.5) to the point where the T'= 1 isotherm reaches x'=l defines A a'. As 
shown in Figs. 3(b) and 3(c), A a' decreases as b increases. In fact, A a' is purely a function 
of b . The dependence of A a' on b , as determined numerically from ( 13), is plotted in Fig. 4 


The cooling incursion A a (in m) is then determined from the inverse of (3e) 


{ TS \ 


A a 


K 


1/2 


33 




x,A a' 


Also plotted in Fig. 4 is an approximate formula for A a' (b ) , given by 


(14) 


A a' 


— [l-tanh(. 116-10)] + 


1.415 

yfb 


[l + tanh(.1046-.4)] 


(15) 


Substitution of this formula (15) and (3c) into (14) then gives 


A a = 332 ( K n K x) _ tanh ( n^_io)]+ 2 

pc„u 


f v V /2 
K 33 x h 


pc„u 


[l + tanh(. 104(? - .4)] (16) 
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where A a (in m) is the cooling incursion inward from the tape edge. For b > 10, the second 
term on the right hand side is dominant, and (16) can be written 

( V /2 

Aa = 4 (b> 10) (17) 

pc p u 

v p 

4. EXAMPLE 

This example illustrates the application of ( 14)-( 16) and Fig. 4. It starts from the present 
operating condition of the Langley ATP robot, where the consolidation length x h = 1.27 x 1 0 3 
m and the speed U = 0.0425 m s' 1 . For PEEK/carbon fiber tape laid-up unidirectionally, the 
thermal conductivity components are given by K n = 6 W m 1 K 1 and K 22 = K 33 = 0.72 W m 1 
K' 1 . The density and specific heat are given by p = 1560 kg m 3 and C p = 1425 J kg' 1 K' 1 . 

Substituting these values in (3c) gives b = 10 . Figure 4 or, alternatively, equation (15) then 
gives A a'= 1. Substitution in (14) or (16) finally gives the cooling incursion inward from the 
edge of the tape as A a = 0.4 mm. This point is plotted on the lower curve of Fig. 5. 

As the performance of the Langley robot is upgraded in the future, both U and x h will 
increase. The increase in x h results from a desire to make the consolidation roller more 
conformable. Also, increasing x h increases the effective bonding time, as show by sensitivity 
studies presented in ref. 10. Suppose that both U and x h are doubled. Then b = 40 and a 
second point can be calculated and plotted on the lower curve of Fig. 5. This entire curve is 
obtained by proportionally increasing (or decreasing) U and x h from the initial point 
(b = 10) . 

If successive layers are angled to make the composite substrate isotropic, the average values of 
the conductivity components become K n = K 33 = 3.36; K 22 = .72 . Recalculation with these 

values (and recalling that b depends on K n ) gives the upper curve in Fig. 5, where the 
minimum cooling incursion Ac/ = 1.2 mm. 

5. DISCUSSION 

A comparison of Figs. 4 and 5 shows that, although A a' is a decreasing function of b , A a 
levels off for b > 10 . Equation (17), which is valid for b > 10, provides the following 
explanation: If U and x h are increased proportionally, as done in the example shown in Fig. 
5, then A a remains constant when b > 10. 

The range (b > 10) is of great interest in ATP. Increases in the speed U benefit the 
productivity of a unit, and increases in x h are a natural consequence of the trend toward using 
conformable rollers. As shown by (17), however, increasing x h has the undesirable effect of 
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increasing the cooling incursion A a . A rearrangement of the plies that increases the average 
value of K 33 also increases A a , as does any decrease in the density p or the specific heat 
C p . To counteract these adverse effects on A a , U could be simultaneously increased, as 
exemplified in Fig. 5. 

Equation (17) shows clearly that for Z? > 10 the cooling incursion A a need not decrease as the 
Peclet number Pe = 2 b = pC p Ux h / K n increases. If the increase in Pe is due solely to a 

decrease in K n , A a will not be changed. If the increase in Pe is due solely to an increase in 
x h , A a will actually increase. Therefore, large values of Pe do not necessarily justify the use 
of a two-dimensional approximation in ATP thermal transport models. 

A simple fix for the tape-edge cooling problem is to extend the thermal input past the tape edges 
by the amount A a , as given by (14), (16) or (17). This approach should help achieve uniform 
strength and crystallinity across the width of the tape. Of course, extending the thermal input 
past the tape edges would cause some substrate reheating as the tape is repetitively placed side- 
by-side. Such reheating must be minimized to reduce thermal degradation. Formulas (14), (16), 
and (17) should help provide these minimum values. 

6. CONCLUSIONS 

The idealized three-dimensional thermal transport problem for the ATP configuration shown in 
Fig. 1 was found to have a relatively simple analytic solution (13) for the dimensionless 
temperature T\x\ /, z') . This solution depends on only two parameters, b = Pe/2 , where Pe 
is the Peclet number, and the dimensionless tape half-width a ' . Plots of the dimensionless 
isotherms on the weld interface reveal premature cooling along the side edge of the tape. This 
cooling results from the large thermal gradient between the newly placed tape and the cool 
substrate beside it. 

The prematurely cooled dimensionless distance inward from the tape edge A a' was determined 
from the T'= 1 isotherm, which corresponds to the glass transition temperature T g . The 

inward sweep of this isotherm at the downstream end of the consolidation head, where x h '= 1 , 
defines A a ' . As determined numerically, A a' was found to be a monotonically decreasing 
function of b alone. 

An approximate formula for A a'(b) was also found (15). This formula was then inverted to 
give the cooling incursion Aa (in m). For the important range b> 10, the formula for A a 
reduces to a simple form (17). This simple form shows, surprisingly, that large values of the 
Peclet number Pe = 2b do not necessarily justify the two-dimensional thermal approximation 
that is often used in ATP models. 

For the Fangley ATP robot, presently operating at b = 10 with PEEK/carbon fiber composite, 
the cooling incursion is in the range (0.4 mm < A a <1.2 mm). The lower value corresponds to a 
unidirectional lay-up, and the higher value to an isotropic lay-up. An upgrade in roller 
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conformability that increases the consolidation length x h would increase A a unless a 
compensating change is made, such as proportionally increasing the material speed. 

A simple fix for the tape-edge cooling problem is to extend the thermal input past the tape edges 
by the amount A a , as given by formulas (14, (16), or (17). This approach should help achieve 
uniform weld strength and crystallinity across the width of the tape. 
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Figure!. Three-dimensional sketch o! the automated tape placement (ATP) process. Heat q Q is 
directed at the nip (the line defined by the z -axis) and distributed uniformly over the width of 
the tape {-a < z < a) . Pressure is applied over the distance (0 < x < x h ) . 



Figure 2. Two-dimensional approximation to the three-dimensional model shown in Fig. 1 
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Figure 3. Isotherms of dimensionless temperature T' on the weld interface ( x'-z plane ) for 
dimensionless half-width a'= 1.5 and various b . The x' -axis (vertical axis) lies along the 
centerline of the tape. The edge of the tape is at z'=a'= 1.5. Consolidation ends at x' = 1 . 



(b) b = 31.6 








